Sound propagation in a Bose-Einstein condensate at finite temperatures 
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We study the propagation of a density wave in a magnetically trapped Bose-Einstein condensate 
at finite temperatures. The thermal cloud is in the hydrodynamic regime and the system is therefore 
described by the two-fluid model. A phase-contrast imaging technique is used to image the cloud of 
atoms and allows us to observe small density excitations. The propagation of the density wave in 
the condensate is used to determine the speed of sound as a function of the temperature. We find 
the speed of sound to be in good agreement with calculations based on the Landau two-fluid model. 
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INTRODUCTION 

Long wavelength excitations of a Bose-Einstein con- 
densate (BEC) with repulsive interactions exhibit a 
phononlike, linear dispersion, causing these excitations 
to move at a finite speed c, the speed of sound. Exci- 
tations with a wavelength comparable to the size of the 
cloud result in collective shape oscillations of the system 
[T]. The dispersion relation is linear up to wave vectors 
of the order of the inverse of the healing length. Pertur- 
bations with a wavelength much shorter than the axial 
size of the cloud and larger than the healing length there- 
fore give rise to the excitation of a hydrodynamic mode 
propagating at the speed of sound. 

In the case of liquid helium sound excitations have been 
extensively studied below the A-point, where a superfluid 
and a normal fluid coexist. In this regime two sound 
modes can be distinguished: the first mode, first sound, 
consists of an in-phase oscillation of the superfluid and 
normal fluid component, while the second mode, second 
sound, consists of an out-of-phase oscillation of the su- 
perfluid and normal fluid component. The occurrence 
of two distinct modes is caused by the presence of both 
a superfluid density and normal fluid density, which are 
coupled. One of the drawbacks of liquid helium is that 
the interactions are so strong that a clear distinction be- 
tween the two components is difficult which complicates 
the interpretation of the phenomena. 

In the dilute gaseous BEC studied here a two-fluid 
system exists below Tc, analogous to the case of liquid 
helium. In our setup the thermal cloud is in the hy- 
drodynamic regime, where collisions between atoms are 
rapid enough to establish a state of dynamic local equilib- 
rium in the non-condensed atoms. In contrast to Bose- 
condensed liquid helium, the superfluid in the gaseous 
BEC corresponds directly to the Bose-condensed atoms 
and the normal component directly to the thermal, non- 
condensed atoms, since the interactions are much weaker. 

First and second sound in a Bose gas exhibit different 
features than those in a Bose liquid j2 . In liquid he- 
lium, the coupling between the density and temperature 



is weak, since Cp ~ C„, with Cp^y the specific heat for 
constant pressure and constant volume, respectively. As 
a result, in liquid helium first sound is mainly a density 
wave, while second sound is an almost pure temperature 
wave. In contrast, in a Bose gas, the density and temper- 
ature fluctuations are strongly coupled, since Cp/Cy 9^ 1. 
The first sound mode in a Bose gas is largely an oscil- 
lation of the density of the thermal cloud (the normal 
fiuid) and second sound is largely an oscillation of the 
density of the condensate (the superfiuid) |3] . In super- 
fiuid helium, only the first sound mode can be excited 
by a density perturbation. In contrast, second sound has 
a significant weight in the density response function in 
superfluid Bose gases at finite temperatures and can be 
excited by a local perturbation of the density. 

Since we can directly image BECs and make a clear dis- 
tinction between both components, it allows for a direct 
comparison with theoretical descriptions of the two-fiuid 
system modeling the interactions between both compo- 
nents. Thus the research using weakly interacting Bose 
gases promises results that will go beyond the results ob- 
tained using liquid helium. 

In a pioneering paper by Lee and Yang the speed of 
first and second sound is derived for a dilute Bose gas, 
although there is no coupling between both sound modes 
[3]. This interaction is taken into account in the two-fiuid 
model developed by Landau for liquid helium and in 
a hydrodynamic model developed by Zaremba, Griffin, 
and Nikuni for trapped Bose gases [Hj. In these papers it 
is shown that the hydrodynamic second sound mode at 
finite temperature extrapolates to the T = Bogoliubov 
phonon mode. 

The propagation of sound in a harmonically trapped, 
almost pure BEC in the coUisionless regime has been ob- 
served experimentally in a pioneering experiment by the 
MIT group [3 [S] and studied theoretically by various 
authors [1 [3 O HUJ [H [12] . After the first experiment, 
sound propagation has been observed for a BEC in an op- 
tical lattice, the excitation spectrum of a BEC has been 
measured and the excitation of shock waves is observed 
[HIIIITS]. 
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The work presented here describes the experimental 
observation of a propagating sound wave in an elongated 
BEC at finite temperatures and extends the study by 
the MIT group in two ways. First, in the work presented 
here the propagation of a sound wave is observed at fi- 
nite temperatures. The thermal cloud is in the hydro- 
dynamic regime above and the cloud is therefore a 
two-fiuid system below Tc: a superfluid BEC coexists 
with a normal fluid of thermal atoms. If the thermal 
cloud is already close to the hydrodynamic regime above 
Tc, it will be deeply in the hydrodynamic regime when 
the BEC forms, since the collision rate 7 is dominated 
by the collisions between the condensed and the thermal 
atoms determined by the rate 712 |16j . In the two- fluid 
system interactions between the superfluid and the nor- 
mal component are expected to play an important role. 
Second, the high signal-to-noise ratio of our imaging tech- 
nique allows us to make smaller excitations than is used 
in the MIT experiments, thereby limiting non- linear ef- 
fects. The large atom number BECs in combination with 
the weak axial conflnement results in typical axial BEC 
lengths of more than 2 mm. This allows for the determi- 
nation of the speed of sound with a high accuracy, since 
the propagation distance can be large before the sound 
wave reaches the edge of the BEC. 

SOUND PROPAGATION IN A DILUTE 
BOSE-CONDENSED GAS 

Since second sound at finite temperatures extrapolates 
to the T = Bogoliubov phonon mode, we start this 
discussion in the T = limit . The speed of second sound 
in the absence of a thermal cloud can be derived using 
the Gross-Pitaevskii equation (CPE). Reformulated as a 
pair of hydrodynamic equations, neglecting the quantum 
pressure and after linearization, the GPE can be written 
in the simplified, hydrodynamic form |12j 

^ = V(c2(r)V5n), (1) 

where the departure of the density from its equilibrium 
density ricq is given by 5n(r, t) = n(r, t) — ncq{r) and the 
local speed of sound c(r) is defined by 

TOc2(r) = M-Kxt(r=), (2) 

where is the chemical potential, Vext is the external 
confinement and m is the mass. In a uniform Bose 
gas, Vcxt — 0, the speed of second sound is given by 
c — yj ^/m — a/ g nc/m, with tic the condensate density. 
This result was first derived by Lee and Yang [3] based 
on theory developed by Bogoliubov [17] and is therefore 
often referred to as the Bogoliubov speed of sound, which 
we will refer to as cb in this paper, cb only depends on 
temperature through the BEC density and is indepen- 
dent of the thermal density. 



The experiments are conducted in an elongated 3D 
trap where the external confinement in the radial (sub- 
script rad) and axial (subscript ax) direction is given by 

ycxt(a;, y, z) = {i4s,dX^ + ^^re.dV'^ + ^L^^) : (3) 

with tiirad S> LOax, and the density depends on the po- 
sition. Although the confinement is highly anisotropic, 
the BEC is not fully in the ID regime, since ^ ^ /i^rad- 
For the radial average density we use the Thomas- Fermi 
(TF) value n{z) — n{0, 0, z) /2, where n{0, 0, z) is the cen- 
tral axial density. The local Bogoliubov speed of sound 
is given in terms of the radial average density by 

cB(z)=y^=y^^P^, (4) 

V m V 2m 

a result confirmed using different methods [H [lOl [HI [12] . 

In our experiment the axial confinement is weak and the 
density varies slowly along the axial direction. 

To our knowledge two theoretical descriptions are 
available in which the effect of the interaction between 
the normal and the superfluid component on the speed 
of sound is taken into account. Zaremba, Griffin, and 
Nikuni have derived the two-fluid hydrodynamic equa- 
tions for weakly interacting Bose gases ^ and use them 
to discuss first and second sound for a uniform Bose gas 
[TH] . We refer to this description as the ZGN model. In 
the same paper, they derive the Landau two-fiuid equa- 
tions for a dilute gas in a complete local equilibrium 
and use these equations to calculate the first and sec- 
ond sound velocities. We refer to this description as the 
Landau model. In Ref. [9 it is shown the ZGN model 
is valid in the limit in which collisions between the con- 
densed and the non-condensed atoms are ignored on the 
timescale of the collective excitation j^/w 0. The Lan- 
dau model developed for dense fiuids such as superfiuid 
helium is valid for dilute Bose gases in the opposite limit 
of complete local equilibrium ^p,/uj — > 00 [9J. Here, 7^ 
is the relaxation rate for chemical potential differences 
between the condensate and the thermal cloud as given 
in Ref. [9 and uj is the excitation frequency. 

We introduce a measure for the hydrodynamicity of 
the thermal cloud in the axial direction 7 = 722/i^ax) 
where the collision rate 722 — rics a Wjoi is the average 
number of collisions in the thermal cloud. Here, the rel- 
ative velocity Wroi = V^'^gxi where — ^J^k^^TjrrvK is 
the thermal velocity at temperature T and m is the mass 
and CT = Stto^ is the isotropic cross-section of two bosons 
with s-wave scattering length a. Furthermore, noff = 
/ n1^{r)dV/ J n{r)dV = n0ox/\/8 for an equilibrium 
distribution in a harmonic potential, where nOcx is the 
peak density. Written in terms of the number of thermal 
atoms A'ox = "-Oox {'^T^k^T/ (m<D^))"^^^ and the geomet- 
ric mean of the angular trap frequencies u)^ = Lu'^^^aji^^, 
this results in 722 = N corner uj'^ / {2'K'^kBT) « 90s~^ for 
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the highest number of atoms and corresponds to a hy- 
drodynamicity of 7 < 10 in the axial direction. We have 
observed the transition from the coUisionless regime to 
the hydrodynamic regime by studying a thermal dipole 
mode above Tc [H]. Furthermore, it is noted in Ref. [TS] 
that the thermal cloud is deeper in the hydrodynamic 
regime when a BEC forms due to collisions with con- 
densed atoms. As a result, the thermal cloud is even 
more hydrodynamic below than it is above the tran- 
sition temperature. 

Both models are used to calculate the speed of first 
and second sound as solutions of an equation of the 
form — Av? + i? = 0, where A and B are coeffi- 
cients which are given for both models in Ref. [9] . These 
coefficients depend on the condensate density n^, the 
non-condensate density ripx and the temperature T. In 
Ref. |9j the differences between both models for constant 
density n = ric + ?t-cx are found to be very small in the 
case of a weakly interacting Bose gas. As a result, the 
transfer of atoms required to equilibrate the condensed 
and non-condensed atoms is found to play a minor role in 
the determination of the speed of first and second sound 

The speed of first and second sound for both models 
is shown in Fig. [T| where close to Tc the first (second) 
sound mode mainly corresponds to the density wave in 
the thermal (condensed) fraction. The figure shows the 
coupling between both sound modes cause the speed of 
second sound to be smaller than the Bogoliubov speed 
of sound. Around T = O.lSTc this coupling results in an 
avoided crossing. Note that the position of the avoided 
crossing in Fig. [T] is outside the region of validity of the 
model, which is valid for k^T ji [3]. 

EXPERIMENT 

The experimental setup used to create large number 
BECs is described in Ref. [50]. In short, ^'^Na atoms are 
cooled and trapped in a dark-spot MOT and transferred 
to a magnetic trap (MT) after being spin-polarized 21J. 
Forced evaporative cooling on these atoms yields roughly 
1 • 10^ thermal atoms around T = T^. In order to prevent 
three-body losses, which limit the density, as well as to 
increase the collision rate with respect to the axial trap 
frequency, we work with axially decompressed traps and 
have reached the hydrodynamic regime in the thermal 
cloud [22]. The resulting elongated, cigar-shaped clouds 
used for the experiments described here have an aspect 

ratio Wrad/t^ax ~ 65. 

The experiments are conducted on clouds at various 
temperatures below Tc. The number of condensed atoms, 
slightly depending on the temperature, is roughly A^c = 
1.7 • 10^ and the BEC density is about 2.5 • lO^^m-^. At 
the lowest temperatures, the BEC has a radial TF radius 
of roughly 22 /im and an axial TF radius of 1.4 mm. 




T/T, 

FIG. 1: (Color online) The speed of first (I) and second (II) 
sound calculated in the Landau model (solid lines) and the 
ZGN model (dashed lines), normalized to cb(0) (Eq. Q) 
as a function of the reduced temperature T/T^. The inset 
shows the same temperature range on a larger vertical scale. 
The densities and temperature used for the calculations for 
0.2 < T/Tc < 0.9 correspond to the intrapolated experimen- 
tal values. For T/Tc < 0.2 ric is kept constant, and riex is 
extrapolated to ncx = for T = 0. 



The clouds are imaged using phase-contrast imaging, 
the details of which are described in Ref [23]. Briefly, 
the atoms are imaged in situ and in contrast to other 
implementations of the phase-contrast imaging technique 
[24 , our implementation does not aim at non-destructive 
imaging, but uses the periodicity of the intensity of the 
phase contrast imaging technique as a function of the 
accumulated phase of a probe beam that propagates 
through a cloud of atoms. Therefore, the intensity sig- 
nal for large enough accumulated phase show rings in the 
intensity profile. The number of rings is a sensitive mea- 
sure for the atomic density and allows us to determine 
the BEC density, the thermal density and the tempera- 
ture within a few percent. The lens setup used results in 
a diffraction limited resolution of roughly 4 /im. 

For the experiments we evaporatively cool the atoms in 
the presence of a blue-detuned focused laser beam aimed 
at the center of the trap, which acts as a repulsive opti- 
cal dipole trap. The beam is focused using a cyhndrical 
lens yielding a sheet of light with a (l/e)-width of the 
intensity of 127 ± 20 /xm. The non-focused direction has 
a (l/e)-width of roughly 5 mm. The light of the dipole 
beam is detuned 20 nm below the ^^Na D2 transition 
and the power adjusted in such a way that the repulsive 
potential has a height of (0.24 ± 0.04) where /i is the 
chemical potential. Due to the large detuning heating 
caused by light scattering is negligible. The laser power 
of the dipole beam is continuously measured and used in 
an electronic feedback circuit controlling the efficiency of 
the acousto-optical modulator (AOM) that deflects laser 
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light into an optical fiber whose output is used for the 
dipole beam. Using this procedure the stabilized inten- 
sity has a RMS fluctuation of 0.1%. Typically, the power 
is adjusted on the order of milliseconds and the power 
is therefore not stabilized instantaneously, mainly due 
to heating of the optical fibers. However, the height of 
the potential is not critical during the first stages of the 
evaporation (tens of seconds), where the thermal energy 
k^T is large with respect to the dipole potential, ensur- 
ing the stability of the height of the potential long before 
the BEC is formed. An alternative way of exciting a 
sound wave is by turning on the repulsive potential after 
the BEC is formed. Since this procedure is immediately 
sensitive to the height of the potential this alternative re- 
sults in less reproducible excitations than the procedure 
used in our experiments. Turning the dipole beam sud- 
denly off ((l/e)-time ~ 250 ns) causes a local dip of the 
BEC density. This perturbation splits up in two waves 
propagating symmetrically outward, both with half the 
amplitude of the initial perturbation. A schematic repre- 
sentation of the excitation procedure is shown in Fig. |2] 




FIG. 2: (Color online) Schematic representation of the exci- 
tation of a sound wave, where the trapping potential, height 
and width of the perturbation are roughly on scale. At t < 
a BEC is formed in the presence of a blue-detuned repulsive 
dipole beam which adjusts the trapping potential. At t — 
the dipole beam is suddenly turned off and the inflicted den- 
sity perturbation causes two density dips to move outward for 
t > 0, propagating at the speed of sound. 

Clouds are imaged at about ten different times after 
the dipole beam is switched off, where for each shot a 
new cloud is prepared, since the imaging scheme used is 
destructive. The initial conditions of the newly prepared 
clouds show only a small variation, since the density in 
the final stage of the evaporative cooling process is lim- 
ited by three-body losses. The density as a function of 
the axial position is determined by making ID fits of the 
radial profile for all axial positions. The fit function used 
is a bimodal distribution which is the sum of a Maxwell- 
Bose distribution describing the thermal cloud and a TF 



distribution describing the BEC. Each propagation time 
for a series is mostly measured twice. The total accumu- 
lated phase as well as the radial width of the BEC can 
be used as a measure for the local density, although the 
width is expected to yield inferior results due to lensing 
effects and the limited resolution in the radial direction. 
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FIG. 3: Density profiles showing the peak density of the BEC 
as a function of the axial position for seven different propa- 
gation times. The unperturbed density profile is subtracted 
from perturbed ones to increase the visibility of the density 
dips. The dotted line is a guide to the eye following the trail- 
ing edge of the dip. The horizontal scale is given in terms of 
the ajdal TF radius Rtt, where Titf ~ 1-4 mm. 



RESULTS & DISCUSSION 

The resulting axial density profiles of the condensate, 
shown in Fig. [3] for various propagation times, clearly 
shows the outward traveling density dips. In the density 
profiles of the thermal cloud these dips are absent, and 
no propagation is observed. The propagation of sound in 
the thermal cloud is discussed at the end of this section 
and for now we consider only the density perturbation of 
the condensed fraction of the cloud. 

The density profiles shown in Fig. |3] show that the 
shape of the perturbation, which is initially approxi- 
mately Gaussian, deforms during the propagation. The 
deformation is caused by the dependence of the speed of 
sound on the density and its effect will be estimated using 
Eq. Q . The estimated change of the sound velocity cor- 
responding to a change in density An is approximated by 
Ac = (c/2)An/n. The dependence of the sound velocity 
on the density causes the center of the dip to move slower 
than the edges. Furthermore, this effects results in the 
trailing edge to 'overtake' the center, while the leading 
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FIG. 4: Simulation of the normalized density profiles An/n 
along the axial axis in terms of the axial TF radius _Rtf for 
different times t after the sound wave is excited. The left 
column corresponds to a weak perturbation An/n — 10~^, 
the right column corresponds to the perturbation applied in 
the experiment, An/n = 0.25 and therefore both columns are 
on a different vertical scale. Note that the t = figures (top 
row) have a vertical scale twice as large as the other figures 
in the same column. The dashed lines in the figures indicate 
An/n = 0. 



edge outruns the center part. The resulting high density 
gradient leads to the formation of shock waves when the 
density gradient is of the order of 1/^, where ^ is the 
healing length [TS]. Formation of shock waves compli- 
cates the propagation due to strong non-linear behavior. 
The typical time for the formation of shock waves is esti- 
mated by considering the difference in traveled distance 
Az of the tailing edge (density n) with respect to the cen- 
ter (density n — An) after t propagation time, yielding 
Az « — (cr/2)An/n [TT]. The edge will reach the center 
when Az ~ a, where a is the width of the perturbation, 
resulting in r > an/{Anc) which is the time for shock 
waves to form. For our typical parameters this yields 
T > 200 ms. In the experiments propagation times are 
less than 110 ms. 

For times less than r we already see that the prop- 
agating wave deforms. We have made a simulation of 
the Gross-Pitaevskii (GPE) equation which describes the 
BEG at T = to analyze the effect the deformation has 
on the propagation of the condensate density wave. The 
GPE is numerically solved using the time-splitting spec- 
tral method described in Ref . [25] . Since the experiments 
are done on very elongated, cigar-shaped BEGs and all 
effects are found in the axial directions, the calculation 
time can be reduced by solving an effective ID equation 
in the limit of strong coupling (26j . Many experiments on 



both the statics and dynamics of BEGs have shown that 
experiments can be modeled accurately by numerically 
solving the GPE, for example in experiments on interfer- 
ometry [27] and superfluidity [55]. Note that we assume 
that the temperature dependence on the deformation of 
the condensate density wave can be neglected. 

The simulated density profiles are shown in Fig. |4] 
where the perturbation is excited by growing the BEG in 
the presence of an additional Gaussian shaped potential. 
At t = this extra potential is suddenly switched off, but 
the harmonic confinement remains. Two situations are 
shown; a small initial perturbation and a perturbation 
corresponding to the situation in the experiments. The 
simulations show the deformation of the initial shape of 
the perturbation and the steepening of the traihng edge 
of the condensate density wave under the experimental 
conditions. After about 200 ms the simulations indeed 
show the formation of shock waves. 

We have run simulations for larger perturbations 
An/n « 0.5 and find shock waves to form within a 
few tens of milliseconds. Gomparing these results to 
the sound propagation experiments by the MIT group 
[71[S], where apphed perturbations are reported as large 
as An/n — 1, suggests shock waves have formed in their 
experiments shortly after the sound wave is excited. The 
MIT group has measured the condensate density wave to 
propagate at cb (0) , even though strong non-linear effects 
are expected to influence the propagation, as already re- 
marked by Kavoulakis et al. |llj . 

In our experiment sound waves are excited in elongated 
cigar-shaped BEGs allowing us to make a perturbation 
much longer than the radial size to ensure the excitation 
of a one-dimensional motion, while it remains small with 
respect to the axial size. In the experiment reported by 
the MIT group the size of the excitation is of the order 
of the radial size of the BEG. 

The simulations for small density perturbations 
{An/n = 10~*) confirm the minimum of the density dip 
moves with cb(0) as given by Eq. ^ and shows the va- 
lidity of the simulation. For the perturbation used in the 
experiment the simulations show the minimum of the dip 
propagates about 8% slower than the speed of sound, in 
agreement with the estimate Ac = —{c/2)An/n. This 
is disadvantageous for the determination of the speed of 
sound from the measured density profiles, since the posi- 
tion of the minimum of the dip is the easiest parameter to 
derive. In order to derive the appropriate speed of sound 
from the density profiles, not only the position of mini- 
mum is determined, but also the amount of deformation 
is taken into account by determining the asymmetry of 
the propagating wave. 

The simulations suggest the density at the (l/e)-height 
of the dips is a good approximation of the unperturbed 
density (c/. Fig. [4|. We therefore use the deformation at 
(l/e)-height in addition to the position of the center of 
the perturbation in these simulations as the measure of 
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the distance the perturbation has traveled. The speed of 
sound is now found to remain constant within 2% for a 
broad range of perturbation amphtudes (10~^-0.3), only 
to show larger deviations when shock waves are formed. 
For typical propagation time used in the experiment, the 
simulations show the decrease of the propagation speed 
due to the axial density dependence remains below 2% 
as well. 

Returning to the results of the experiment, the mea- 
sured density profiles are fitted to an asymmetric func- 
tion, yielding both the deformation and position of the 
propagating condensate density wave. Taking the same 
combined measure for the traveled distance as used in 
the simulations, the propagation distance varies linear 
with the propagation time, the slope of which is used to 
determine the speed of sound. 

This procedure is repeated for different temperatures, 
which are created by adjusting the final rf-field frequency. 
Each series, consisting of about ten shots, is used to de- 
termine the propagation speed for that temperature. Rf- 
induced evaporative cooling does not allow to cool to 
temperatures below k^T w fi, since around this temper- 
ature both thermal atoms and condensed atoms are re- 
moved from the trap by the rf-field, setting a lower limit 
to the temperature reached in the experiment. The high- 
est temperature used in the measurements corresponds to 
T/Tc = 0.72 ± 0.04. For higher temperatures the axial 
length of the BEC becomes more than 15 % shorter than 
the typical BEC length. By choosing this upper limit 
for the temperature prevents the need to account for the 
axial density dependence. 

An unperturbed cloud in each series is used to derive 
the temperature, thermal density and BEC density for 
that series. In these measurements the chemical potential 
based on the axial size of the BEC and the total accu- 
mulated phase agree within 5%. The measured speeds 
are normalized to cb(0) based on the measured BEC 
density. The resulting normalized speed of sound as a 
function of reduced temperature and thermal density is 
shown in Fig. [5] (a). From these measurements we con- 
clude that, even though the thermal density is varied over 
more than an order of magnitude, the effect of the ther- 
mal cloud appears to be constant within the accuracy of 
the measurement. Furthermore, the measured speeds are 
found to be approximately 7% lower than cb(0) given by 
Eq. 0. Note that since the condensate density is limited 
by three-body losses, the variation of function of 

the temperature is modest. 

In Fig. |5] (b) the measured speeds are normalized us- 
ing the speed of second sound calculated using the Lan- 
dau and ZGN model. We find the measured propagation 
speed to be in excellent agreement with the speed of sec- 
ond sound given by both the Landau and ZGN model 
within the accuracy of the measurements. This result 
shows we have measured the effect of the thermal cloud 
on the propagation of a sound wave in the BEC. How- 
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FIG. 5: The normalized speed of sound as a function of the 
thermal density. The upper axis gives the reduced tempera- 
ture r/Tc for the corresponding data point. Figure (a) shows 
the speed of sound normalized to cb(0) (Eq. Q) based on the 
measured BEC density. Figure (b) shows the same measure- 
ments normalized to the speed of sound based on the Landau 
model (diamonds) and ZGN model (squares) based on the 
measured BEC and thermal densities. The main contribution 
to the error bars is the uncertainty in the measured densities. 



ever, the dependence of the speed of second sound on 
the temperature is modest in the experimentally accessi- 
ble temperature range. Furthermore, since the difference 
between the Landau and the ZGN model is smaller than 
the experimental uncertainty we cannot distinguish be- 
tween both models in the current experiment. 

We have not been able to observe a sound wave in 
the thermal cloud (first sound) in this experiment. For 
the higher temperatures this can be explained by the 
moderate perturbation depth with respect to the thermal 
energy. At the lowest temperatures, still above k^T = 
II, the thermal density is small and the thermal cloud 
spatially barely extends further than the BEC, causing 
the signal-to-noise to be insufficient to see small density 
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perturbations in the thermal cloud in this regime. Above 
Tc we are able to observe a density perturbation when 
the excitation is of the order kuT, but the excited wave 
damps too fast to be able to observe its propagation. 

CONCLUSION 

In conclusion, we have excited a hydrodynamic mode 
in a BEC: a propagating sound wave. We measure its 
propagation speed, which is used to determine the speed 
of sound in the BEC as a function of the temperature. 
The combination of the phase-contrast imaging technique 
and elongated large atom number BECs allows us to 
make a moderate density excitation in the BEC and ob- 
serve the propagation of a sound wave. Numerical sim- 
ulations are conducted to model the non-linear propaga- 
tion. We find the speed of sound to be in good agreement 
with both the Landau model and ZGN model in which 
the coupling between the first and second sound modes is 
incorporated and thus we have observed the effect of the 
thermal cloud on the speed of sound in the BEC. How- 
ever, the effect does not vary within the accuracy of the 
measurements with the temperature for our experimen- 
tal conditions. The Bogoliubov speed of sound, in which 
this coupling between the two sound modes is absent, de- 
viates significantly from the measured speeds reported in 
this paper. 
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